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IMPURITY-ACTIVATED CRYSTALLINE PHOSPHORS 
THEIR PRODUCTION AND THERMOLUMINESCENCE CURVES 


By А. WRZESIŃSKA 


Institute of Experimental Physics, Nicholas Copernicus University, Torun 
м 
(Received November 22, 1955) 


Some useful observations made when searching for methods of making strongly 
glowing Phosphors (of not very pure materials) are reported. A simple apparatus has 
been constructed which traces the glow curve directly on the photographic paper by 
means of a light spot. A uniform increase of temperature is easy to obtain. Some ob- 
servations concerning the influence of Mg on the depth electron traps have been made. | 
Two new methods of finding the depth of electron traps (without approximate calcu- 
lations) have been applied to the produced phosphors. 


I. Introduction 


This paper concerns those crystalline phosphors for which a necessary condition 
of glowing is an admixture of impurities, called activators. The activators are usually 
some heavy metals. Natural impurity-activated crystalline phosphors have been known 
ever since the 17-th century, while methods of their artificials production have been 
studied thoroughly since the 19-th century. In this. respect the greatest credit goes 
to Lenard and his school. Later on several investigators worked out different methods 
for producing these phosphors. This problem is still of interest on account of the 
great variety of phosphors and their increasing practical importance. The investigation 
of their phosphorescence and photoconductivity is also important for the understanding 
of their structure and the mechanism of photoluminescence. In particular, the investi- 
gation of glow curves enables the depth of electron traps to be determined, i. e., the 
difference of energy between the localized level and the bottom of the conductivity 
band. The mean time interval during which electrons remain on these levels determines 
the mean М of phosphorescence and іп connected with the depth of electron 
traps by — | 

1 E 


(1) | j om ager 


where тіз the. mean life-time, pg is. а constant which is characteristic for а given 
phosphor, and e is the depth for the trap. Numerous investigations have shown that 
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some of these traps are activator-independent traps and are connected with natural 
defects of the cristalline structure, while the remaining activator-dependent traps are 
due to the activator, or to the so called co-activator. 


II. Production of impurity-activated cristalline phosphors 


The first part of this paper is devoted to methods of producing cristalline phosphors 
with long life-times. Such phosphors are best suited for investigations of glow curves. 
The phosphors produced from: CaS, SrS, and BaS have particularly lasting phospho- 
rescence. The methods of their production are well known and described in the lite- 
rature, but the utmost purity of materials is necessary for their production. In some 
cases even a purity above the “spectrally pure“ is needed, otherwise the intensity of 
phosphorescence is considerably reduced. Such substances are not accessible and 
suitable quartz and platinum apparatus for purification was not available. Neverthelless, 
we attempted to produce phosphors from the available materials viz. from СаСО;, 
SrCO,, and ВаСО,, sulphur, several salts of light metals for melting, as well as the 
metals Bi, Cu, Tl, Mn, Zn, Pb as activators. 

As a rule we followed the methods to be found in the literature especially those 
of Lenard and Tomaschek; but searching for strongly luminescent phosphors, derivable 
from the available substances, we noticed some useful facts: 

(i) Some of the melting materials recommended in the literature may be replaced 
successfully by some others. For example, when producing the BaS-Cu phosphors, 
we used Na,B,O, for the melting instead of the recommended Li,PO,, and obtained 
a more strongly glowing phosphor. 

(ii) While making several trials, we found that the way in which the firing of the 
phosphor is carried out is of importance for the strength of luminescence. In particular, 
if the substance to be heated is put into the oven at a comparatively low temperature 
(е. в. CaS—Bi at 400°C) so that the increase of temperature up to about 1100°C is 
slow, we get a weakly phosphorescent-low quality phosphor. 

(ii) An admixture of MgCO, and an admixture of glucose to the mixture to be 
heated usually increases the intensity of luminescence. 

For heating of phosphors, quartz crucibles, or, more frequently, crucibles of 
glaze porcelain were used together with a suitable layer, to prevent contamination of 
_ the phosphor by the glaze. After several trials, we obtained the following phosphors 
which are distinguished by intense luminescence. 

1) CaS—Bi luminescence violet 

2) The mixture CaS — BaS — SrS — Bi, green or blue зоа to the proportion 
of the three basic substances. 

3) BaS—Cu orange 

4) Са5--ТІ green, weaker than the above mentioned phosphors. 


These phosphors glow strongly when excited by the light of a mercury lamp; 
‚ but phosphors 1 and 2 are also good if excited hy daylight: 
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Ш. Glow curves 


A. The second part of this paper concerns the investigation of the glow curves 
of the produced phosphors. An estimate of the depth of electron traps is of theoretical 
value, since it may elucidate the character of the structure defects causing the appear- 
ance of those levels. 

Moreover it is of practical value for applications involving phosphors to know the 
position of the glow-curve (thermoluminescence curve) maximum. 

The glow curve is the plot of intensity of luminescence versus temperature, where 
the temperature is increasing linearly in time. In order to obtain entire glow curve 
it is necessary to excite the phosphor at a sufficiently low temperature, so that the 
phosphorescence is completely frozen-in. Then, after interrupting the excitation, one 
measures the change of intensity by increasing the temperature linearly in time. 

B. Apparatus. Our apparatus was constructed so as to enable the temperature 
to be increased continuously and uniformly, and to allow simultaneous measurements 
of temperature and relative intensity of luminescence. The arrangement for changing 
temperature was similar to that used by Randal et al. (1945). It consists of copper 
box with a channel for liquid air surrounded by an electric heater coil. Both are thour- 
oughly isolated from the outer walls of the box, except for one of the vertical walls, 
which is in contact with liquid air and with the heater сой. The phosphor was attached 
to this wall (Fig. 1). 


Several ways of fastening were tried (e. g. by means of water glass or polystyrene 
in benzene. No significant influence of the glue was noticed. Usually however, the 
phosphor was held on the wall without glue, with the help of a quartz window and hard 
strings. One of the thermocouple junctions was fastened to the same wall. It should be 
placed so as to indicate the actual temperature of the phosphor without making it 
difficult to press the phosphor to the wall. The junction was flat and situated in a shallow 
cavity in the middle of the wall. The current produced by thermocouple passed through 
a properly shunted galvanometer. The shunt was adjusted so that the change of tempe- 
rature AT = 1° corresponded to a displacement of the light spot (temperature spot) 
by 1 mm on the scale. The window pressing against the phosphor froze during excit- 
ation. It is evident, that the most practical means to prevent this is to use a double 
window in form of a cylinder 25 mm long with a cavity into which small quantity of 
Р.О; was placed to remove the moisture. 


The apparatus for the measuring the relative intensity of phosphorescence consist- · 
ed of a photocell, a d.c. amplifier and a mirror galvanometer. The phosphor was cooled 
up to about 130°K and then excited with a mercury arc. After switching off the lamp, 
the phosphor was turned towards the photocell by a special arrangement. The shutter 
of the metallic shield to the photocell was opened for 5 sec. after the excitation was 
switched off. The amplified current was detected by the mirror galvanometer. The 
light spot from this galvanometer was reflected on to a screen placed at a distance of 1 m 
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and indicated the relative intensity of phosphorescence. The remaining part of the appa- 
ratus was constructed so that the glow curve was traced by the light spot directly on the 
photographic paper moving uniformly with the screen in a vertical direction. At the 
same time, the temperature spot traces the change of temperature in time on the same 
paper. The photographic paper was placed under a transparent scale. In order to obtain 
a true glow curve, the temperature of the phosphor must increase linearly with time. 
In order to secure this, a curve (corresponding to the curve gauging the thermocouple) 
was drawn on the above mentioned transparent scale. In order to obtain a uniformly 
rising temperature, it was necessary to adjust the heating current, thus constraining 
the temperature spot to move along the curve drawn on the scale. Obviously, this curve 
has to be adjusted to the desired speed of heating, whereas the speed of the screen is 
fixed. The speeds of heating used were usually 1 grad/sec. and 0.5 grad/sec., while 
the speed of the screen was 0.5 mm/sec. The temperature range for which the glow 
curves were obtained was 130°—600°K. After the two spots have traced the intensity 
and temperature curves, the entire screen was exposed to the light for a few seconds, 
in order to obtain the negatives of the scale and the temperature pattern curve as well. 

The rate of temperature rise was regulated by means of a resistance in the heating 
circuit. The role of the experimenter is to observe the temperature spot and to regulate 
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Fig. 1. Block diagram of the apparatus. Hg — mercury lamp, O 1 i 

. c i : y p, О — electric heater coil, T — thermo- 

couple, Р; — investigated phosphor, Р, — photocell, Р, — phosphorus-pentoxide, E C. 4. — D.C. 

Ampifier, Сі, С, — mirror galvanometers, К, — intensity light spot, R, — thermocouple light spot, 
251, 5, — bulbs. К — opening of the channel for liquid air, L — axis of rotation. 
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the resistance so that the spot moves along the pattern curve. After a certain amount 
of training this could be achieved quite satisfactorily, as may be seen on the curves of 
Fig. 2 and Fig. 3. The whole apparatus was operated by a single person. 


Fig. 2. Example of a glow curve for CaS-Bi (with ^ Fig. 3. Example of a glow curve for green phosp- 
Me) ее range 180° К —570^K. Original. hor. Temperature range 135° K—570°K. Original. 


600 500 400 30 200 {00 | 
Fig. 4. Glow curve of violet CaS-Bi phosphor (without Мр); Temperature range 130°К—600°К. 


С. Discussion of glow curves. The glow curves of two kinds of Саб — В! 
phosphors have been investigated. One of them was produced by heating at a tempe- 
rature of 1100°C the mixture: 6g CaCO + lg S + 0.25 glucose T 0.2g Ха,50, ay 
0.2g K,SO, + 0.0002g Bi, the other was produced by heating under the same condi- 
tions a mixture which contained additionally 15 MgO. It appeared that the mixture 
containing this last substance increases the brightness of the phosphor and causes 
considerable changes in the distribution of the trap depths. The glow curve of the so 
called *normal* CaS—Bi phosphor may be found in the literature (Garlick 1945). 
It appears that it is precisely the thermoluminescence curve of our phosphor obtained 
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without the mixture of MgO. This curve is given in Fig. 4 for comparison with the ther- 
moluminescence curve for the same phosphor, with the mixture of MgO (Fig. 5). 
The curve (4) like the curve in Garlick’s book for the “normal CaS — Bi phosphor 
possesses three maxima. The first of them, a very weak one, corresponds to the tempe- 
rature 153°K. The second about 390°—400°K is the strongest maximum and is 


J 


600 500 400 300 200100 


Fig. 5. Glow curve of violet CaS-Bi phosphor (with Mg). Temperature range 130° К — 600° К. 


M 


600 500 40 %0 200 00 
Fig. 6. Glow curve of green CaS-SrS-BaS-Bi phosphor. Temperature range 130° К — 600^ K. 


rather diffuse. The third one, about 553°K, is nearly as strong as the second. Fig. 5 
represents the glow curve of the CaS — Bi phosphor with admixture of Mg differing 


widely from the previous case. Besides a weak maximum at about 153°К there exists 


also a new weak maximum at about 258°K. The third maximum about 390° —406°K 
is much stronger than the corresponding maximum of Fig. 4. On the other hand, 
the maximum at the high temperature (of 553°K) is missing. It is seen that the admixture 
of MgO considerably influences the structure of the electrons traps, i.e., their depth 
(appearance of new maxima and disappearance of other maxima) as well as the concen- 
tration of traps of a given depth (the increase of the intensity at about 400°K). 

Since the maxima at about 255° —270* К appear also for other phosphors with 
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an admixture of Mg (Fig. 6 and Fig. 7) it seems that Mg plays an essential role in the 
forming of traps corresponding to this maximum. 

The glow of green and blue phosphors are very similar and possess three maxima. 
The main difference is that in the former, the maximum appearing at the highest 


500 400 300 200 100 
Fig. 7. Glow curve of blue CaS-SrS-BaS-Bi phosphors. 


600 500 400 300 200 100 


Fig. 8. Glow curve of orange BaS-Cu phosphor. Temperature range 100°K—500°K. 


temperature is shifted by about +10°K. Perhaps it is connected with the fact that 
the green phosphor contains more CaS. This supposition is plausible, since the violet 
phosphor BaS —Bi (whose basic substance is CaS) possesses a maximum ata still high- 
er temperature. 

The BaS—Cu phosphor orange, possesses no high temperature peak, Therefore 
its life-time at normal temperature is shorter than the case of the above mentioned 
phosphors. Its maxima correspond to the temperatures 150°K and 296°K (Fig. 8). 

D. The depth of electron traps. The depths of traps have been determined 
from the temperatures of the maxima of the thermoluminescence curves. To this 
end was found the extremum of the theoretical curve derived under the assumption that 
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the retrapping effects (Randall at al. 1945) does not exist. Although the neglect of 
retrapping does not seem to be well founded (at least in the case of photo-conductive 
phosphors) nevertheless, it appears to have little influence upon the temperatures at 
which the peaks of glow curves appear. The error caused by the neglect of retrapping 
in the computation of е does not exceed 3 per cent i. e., it lies within the limits of 
experimental error. Under the assumption that all traps have the same depth the 
expression for the glow curve without retrapping 


ғ € 


Тр 
(2) Тыс МЇ, 86 RT ре КЕ 
where C fixes the unit fot Г and may be taken С = 1, T, is the temperature of excita- 
tion, ло is the density of the excited centres at the moment the excitation is switched 
off, В Һе speed of heating phosphor in grad/sec, T—the temperature at a given 
instant. The condition for the maximum of the curve (2) yields 


= en KT, 
В 


from which 


where Т, is the temperature of maximum of the curve. Thus, 


x fue De 
(3) кг tex n KT 
` From (3), е may be determined by a graphical method. Пен 
RR OUT M — Ро Ты 
ЕТ ЫТ 


we have 


ge =U ora ME INU 


From the curve z + іп z versus z, we find z and e for a given U. In order to know U, we 
must know ро for a given phosphor. According to Randall, p, is of the order 108*1 зек-1 
for phosphors whose basic materials are alkaline earth sulfide. In this paper the value 
Ро = 10% зеК 1 has been assumed. We have also applied another method for computing 
the depth of traps, In this method, it is necessary to know only two temperatures, 7) 
and T, of the same maximum corresponding to two different speeds of heating 8, and 
Вә. Since the depth is independent of the heating speed, we find from (3) 


КУЕП Ке: 
ТЫРЫ т Ту ae 
4 € б 3 е В, 
whence 
Т? 
1; des E 7, ln — 

(4) | In РК _ A | B. 

. A 46 Ti өз ipe 
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Denoting 
K 
In ES 
& 
we have 
Е. = 
(5) Po = K e? 


Again, substituting pọ from (5) into (3) we find for = 


T 2 
6 a OF In ez a 
(6) KT, (1 в + | 


The idea of finding py from the difference of position of the maximum at two different 
speeds of heating is due to (Parfianowicz 1954). However, the author has given for the 
computation of the depth of traps only an approximate formula, whereas, by a sub- 
sequent application of this method, we find an expression for є free of any approxima- 
tions. 

The depths of traps coniputed by the two methods are given in Table 1.. The 
results agree with 12 per cent. 

E. Discussion. In the first method, the error of the value of & is due mainly 
to the circumstance that the constant ру is known only up to the order of magnitude 
po = 1081зеК-1. The error of e does not exceed 12 per cent. On the other hand, the 
experimental error in the determination of the temperature of the maximum 47, = 
= 2°К. is negligeable here and causes Ae — 1 per cent. 


Table I 


В, == 0.5 
grad/sec 


В, =1 


" Method I py = 108 вес"! 
grad/sec 


0.269 eV 
0.465 eV 
0.747 eV 


0.269 eV 
0.707 eV 
1.103 eV 


CaS —Bi 
without Mg 
(violet) 


CaS— SrS — BaS— Bi 0.278 eV 0.299 eV 
(green) ; 0.488 eV 0.55 eV 
0:677 eV 0.60 eV 


0.273 eV 
0.478 eV 
0.659 eV 


0.258 eV 
0.537 eV 
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In the second method the main error for 8. is connected with the error in deter- 
mining the temperatures of the maximum T, and Т,. The quantity a, to be computed 
from (4), is very sensitive to these errors. For diffuse, high-temperature peaks, this error 


ay В ; а 
amounts to 4Т--29К in the conditions of our experiments. As e =КТ, | In — + a 


a simple computation yields, in the case ДТ, = 2°K. for T in the interval 350° —550°K, 
again an error for є of 20 per cent. 

At low temperatures the error for e. would be greater for the same AT, : however, 
in this case the peaks are not so broad, and AT, = 1°K. Thus, in both cases, the values е 
are subject to an error of about 20 per cent. 

Let us discuss now the errors due to the neglect of the retrapping of electrons. 
Assuming that the retrapping actually takes place (kinetics of II-nd kind — Curie 1955) 
and that the probabilities of the retrapping, and return of the electron to the centre 
of luminescence, are equal, one finds a formula 


Dee Re. 
(7) J= e pa 
— 0 - 5 
ERES 


(see Romanowski 1946) 


where v is the density of traps. The condition for the maximum (by exciting the phos- 
phor to saturation — Parfianowicz 1954) yields 


2 
(8) АКТ т 2 ро Tn? К 
ep 
whereas, by neglecting the retrapping it was 
a a 
б == KUL, pie 
whence 
2 
(9) Da KEE (ele + № 2) 
ep 
5 Ро f Met K 2 5 
Since In E dis of the order 20 and іп 2 of the order 0.7, е” is larger than ғ by 
3.5 per cent. 


On the other hand, the value £ computed by the second method is independent 
of whether one starts with formula (2) or with formula (7) (corresponding to the 


kineties of the first or second kind). Applying (8) to the computation of a’ = In К | 
Е 
ме һауе | 
T T2 
р м 
(10) a’ Жы Coh 737241 


„Тү Т; i 
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Comparing (10) with (4) we find a’ =a—In2 


but 
T 
FEE RT Раса 
p E 
or 
, 26 ГД 
(11) E = +a + In2 


Comparing (11) with (6) we see that є' = e 
Accknowledgement. The author is indebted to Professor A. Jablonski for 
suggesting this investigation and for his interest in this paper. 


KPATKOE СОДЕРЖАНИЕ 


A. Вржесинская, Активированные кристаллические фосфоры; ux произвотдство 
и кривые термолюминесценции. 


В поисках методов получения ярко светящихся фосфоров из находящихся в продаже 
материалов произведены следующие наблюдения: а) Чтобы получить ярко светящийся 
фосфор, следует материал, предназначенный для обжигания, вложить в печь, имеющую 
уже заранее соответствующую температуру, чтобы получилось стремительное, а не мед- 
ленное, его нагревание. в) Некоторые плавни, имеющиеся в рецептах в литературе, мож- 
HO с большим успехом заменить другими. Получены следующие ярко светящиеся фосфо- 
ры: г) Са5-В:, светящийся фиолетовым светом, 2) фосфорная смесь в составе: CaS, SrS, 
BaS-Bi, светящаяся зелёным светом или голубым, в зависимости от взаимоотношения 
выступающих в ней основных вещесть, 3) фосфор Ва$-С:, светящийся оранжевым светом, 
4) фосфор Са5-ТІ, светящийся зелёным, более слабым, чем предыдущий, светом. 

Для исследования кривых термолюминесценции сооружена такая аппаратура, бла- 
годаря которой кривая накаления чертится сразу посредством светильно-напряжённого 
пятнышка на фотографической бумаге. Однообразие возростания температуры контроли- 
руется специальным устройством. При помощи этой аппаратуры получены кривые 
накаления фосфоров, изготовленных в первой части этой работы. Произведены некоторые 
наблюдения относительно роли примеси магния в структуре электронных ловушек. 

Глубину ловушек находили с положения (в температурном масштабе) максимум кривых 
термолюминесценции. Применены при этом два метода. При первом находими глубину 
ловушек, решая уравнение на максимум кривой термолюминесценции способом графичес- 
ким. (Нужная для конкретного исчисления глубины ловушки величина постоянной, харак- 
терной для определённой группы фосфоров, взята из литераруры. Здесь дело в постоянной 


Ро› выступающей в формуле для вероятности освобождения электрона из ловушки по- 
€ 


средством тепловой энергии: E. p. = poe- KT 
Метод второй состоит в применении (внушенного Парфиановичем 1954) метода отыска- 
ния постоянной p, из двух температур T, и Ta, отвечающих тому же самому макси- 
мум кривой термолюминесценции при разных скоростях нагревания. Применяя последо- 
тельно этот метод, получена точная формула для исчисления глубины ловушки. Резуль- 
таты, полученные посредством обоих методов сходны в пределах экспериментальных 
` ошибок. Оказывается при этом, что величины для глубины ловушки, полученные мето- 
дом впервым, не зависят от того, будет ли взята для кривой накаления формула, отвеча- 
ющая кинетике первого вида (кинетике мономолекулярной), или кинетике второго вида 
(кинетике бимолекулярной). 
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ON RELATIVISTIC WAVE EQUATIONS WITH A MASS SPECTRUM! 
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The question of relativistic wave equations having a mass spectrum and containing 
new continuous variables is discussed. Mostly the case for which the new varia- 
bles are components of some space-like four-vector is considered, since similar equa- 
tions have been very often used in recent times in connection with attempts to build 
a non-local field theory. 


In recent years, in a number of papers on non-local field theory, relativistic wave 
equations with a mass spectrum and containing new continuous variables have been 
discussed. Usually, however, no attention is paid to a number of difficulties and 
characteristic peculiarities connected with the use of these equations. This is why it 
seems to us worthwhile to publish the present paper devoted to equations of the 
above-mentioned type. 

After the introductory § 1, we will examine in § 2 a number of relativistic wave 
equations with a mass spectrum which have been used till now. In § 3, some of their 
general properties will be discussed. 

§ 1. Already about 15 years ago, on the basis of some physical considerations 
which we will not discuss here, (see [1]), investigation was begun of relativistic wave equa- 
tions with a mass spectrum. The latter means that the equations have solutions corres- 
ponding not to one, but to two or more values of rest mass. The equations considered 
at first were G- 5) апа (1-2), ротор to a particle which can be found in 

. two states with pest masses and „with spins = 1 and 3 ; respectively, or 1 and 2 respec- 
tively. The first of these equations is of NUES to ses as a generalized model of a nuc- 
leon and was therefore more extensively investigated [1, 2, 3]. Subsequently, equa- 

. tions in general form of the type 


әт 
Pies. ; = 1 
АИ 0 (1) 


1 This paper has been prepared оп the basis of a lecture delivered by the author on April 4, 1955, 
at the All-Union Conference on Quantum Electrodynamics and Theory of Elementary Particles. The 
appearance of this paper in a Polish Physical journal seemed appropriate, because of the great attention 
paid by Polish physicists to non-local field theory and to the above-mentioned class of relativistic wave 
equations. 
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were studied [4], where L, is a finite or infinite matrix and x, а number different from 
zero. Regarded as equations of this type САС the widely used equations for particles of 
spin 0, 1, equations for particles of spin 1, 5 3 and higher (5, 6], as well as equations for 
particles with two values of spin. The р investigation of equations of type (1) 

revealed a number of interesting features, but in the case of a large or infinite number 
of values of rest mass, the use of such equations (or more accurately, systems of a large 
or infinite number of equations) for any practical computations is not possible. In 
additions, no less important is the fact that in (1) the choice of matrix L,, which 
determines the values of the rest mass, is left most arbitrary. Therefore, it was natural 
to seek a constructive path for finding the relativistic wave equations with mass 
spectra without the use of matrices, which are very cumbersome to work with even 
for the equation (1—3). Such a non-matrix approach was already outlined at the end 
of Reference [1] and as a concrete example was developed in [7] in connection with 
the introduction of new independent variables in the wave function. In other words, 
the function y is considered to be dependent not only on the coordinates and time ж; 
(i = 1,2, 3, 4) and spinor indices, but also on some other four-dimensional quantities, 
for example, some vector и,. It is immediately obvious that introduction of similar 
continuous variables is equivalent to the consideration of infinite-dimensional quanti- 
ties such as р in (1) when Г, is an infinite matrix. A function of the type У (x; ш) 
can, of course, already be considered a scalar, or in the extreme case, can be consi- 
dered as having a small number of spinor indices (bispinor, vector), since the simul- 
taneous use of both vector и; and larger (infinite) numbers of indices would be con- 
trary to the very purpose of introducing now continuous variables. 

In [7], as we shall see, nearly all the typical peculiarities of the equations of the 
class under consideration were elucidated, and by virtue of the difficulties encoun- 
tered, we practically stopped further work in this direction. However, in a whole 
series of papers [8— 14] which have subsequently appeared, discussion of the relati- 
vistic wave equations for the function W (a;,u;) began а new, but usually without 
taking into consideration the difficulties we have referred to. On the other hand, one 
of the possibilities disregarded in [7] is, perhaps, of interest, and does not deserve 
to be neglected. All this prompted the author to return once again to the discussion 
of relativistic wave equations with mass spectra. 

8 2. Let us assume that the wave function V/ is a scalar dependent on the variables x, 
and u;, that the wave equation is not higher than of the second order with respect to 
x;, and finally, that the wave equation does not, of course, contain x; explicitly. Then 
this equation will have the form 


8% 9 9? o9 8 
= Фф i Uili ^, =, — қалын» = 
(за pe, a E Ы дш Ju; ди; ' Эм Эщ eg pres „= 2 


where Ф is some function of its arguments which satisfies the requirements indicated 


above. (In some cases supplementary conditions of the type (2) can be imposed on 
the function). 


| 
| 
: 
| 


On relativistic Wave Equations with Mass Spectrum 165 


In [7], only equations of the type (2) not containing the quantities ши; were 
considered. In connection with these the operators 


д д д 
Da Nc = ig deer 
к= өш > و‎ eae (Mass) @) 
Mig Moog Mu Mig = Ma Sim + Мы ды Е Min ды + Mont д, (3a) 
were used. 
We introduce the coordinates 
uj =r ch y sin 9 сов p, из =r ch y sin sin ф, из = r ch y cos Ó 
u, =i rsh y, — оо < y «oo, 0 < «mn, 0 <фҳ< 2л (4) 
where the vector u; is considered space-like (for the possibility of choosing u; as 
a time-like vector, see discussion below). 
In coordinates (4) 


а а 
б ал oats EEE 
duci z E E (sm ? 5) 2: 
аши r a Mi Ми = وو ہس ہس‎ 
Риш За = (5) 


It is obvious from (5) that the operator  M;, M,, does not contain the variable г, 
by virtue of which equations containing this operator are known to be compatible 
with the supplementary condition (и; и; — г?) ¥ (х; и). The same is true of the 
operators M,, themselves, so that in the sequel in connection with equations contai- 
ning u; and 9/9 и; only in the form М, the above mentioned supplementary condi- 
tion will not even be explicitly mentioned. 

First of all, we shall examine the equation 


k Ам, M, xs | 
Еа ео Ma ма) ¥ (su) =0 (6) 


In this сазе the variables may Бе separated and obviously 


У (х, ш) = v (x) 9 (ш), n IA n) V, (xj) = 0 (7) 


1 
Lip = = Ма Mr o = Аф (8) 
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and the rest energy Ey == m, is equal to (everywhere below h/2 л = c = 1) 
Еу? = то? = х? — 410 (9) 


As solutions of equation (8) we look for functions whose square is integrable on 
the hypersurface г? = const, i. е. we require that the integral f; |p|? dS be finite, where 
dS = ch? y sin 8 dy 49 dg is an element of area оп the hypersurface. In [7] such 
a solution of equation (8) was found directly, but we have no cause to dwell on this, 
since the eigenvalues of equation (8) can be found in still more general form by using 
the results of Reference [4]. As already shown in [7], the eigenvalues A, and A, of the 
operators L, and Ly = — 1i е,ь,„ М, Mim» Where г, is the completely anti-symmetric 
unit pseudo-tensor, are equal to 


a = — (2-10) +95 
i= са. (10) 


where j is an integer or half-integer and а is a real, positive number, and where for 
every j there are an infinite number of solutions differing by whole numbers / 2 j, 
|т| < IL (here 1 (I + 1) and m are the eigenvalues of the operators — (Mog + Мал + 
+ Мы? and i M,,). It is assumed that the operators М,, are in reality operators sat- 
isfying relations (3a) and that they can be obtained in various ways. In the particular 
case where they have been obtained! in the form (3), it is easy to see that L = 0 and 
therefore for j 40, a = 0 and that a 52 0 only when / = 0. Further, in the case of the 
space-like и;, the value J = 0 is not permissible (the square of function ¢ is not in- 
tegrable) and thus the solution of equation (8) is 


а E 1 - a 


On the other hand, if vector и; is time-like, equation (8) will have only the solution 
J = 1 = 0, which corresponds to a continuous spectrum A, = a. Below, we will 
always assume that vector u; is space-like. 

Returning to equation (6), we see that the mass spectrum (9) has the form 
Ев? = т? =x? + В (]2—1) апа is infinitely degenerate (since mg does not depend оп 
the number 7 œ j). Solutions of equation (8) are connected with infinite dimensional 
unitary representations of the Lorentz group which are considered in detail in Ref. [4]. 
Since the number l is associated with the spin of the particle, the degeneracy 
of the J—spectrum means that the spin of the particle in a state with a given mass. 
is indeterminate and that equation (6) obviously cannot provide a true picture. There 
а > 

_1 Аз example of а realization different from (3) the following опе may be mentioned 
Ma = Su z ==) Sey = 


where S;, is an antisymmetric tensor. 
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are two ways in which the degeneracy can be removed. The first consists of imposing 
upon У” some supplementary condition such as the equation used in Ref. [7]: 


| 92 | 9 | 
ІШ Мы DIM n (o + 1) Dx x У (xi, ш) = 0 (12) 


In the rest system, this equation has the form [My My — (jp + 01У = 0 and its 
solution satisfying (6) is different from zero if 


+1 + j — =0, 1S j,j 40; land j — integers. (13) 


Equation (13), in turn, has a solution only when jọ which was introduced in (12), 
is an integer; in this case, for a given jọ there exists one or more solutions corres- 
ponding to different spins and rest masses (for details, see Ref. [7]. The rather chance 
character of the spectrum obtained and, above all, the difficulties of introducing 
the interaction — of which we will speak in 8 3 — make it of greater interest to inves- 
tigate the second path. This н consists of introducing into the equation mixed 


terms containing M;, and ——.Thus in [7], there was a discussion of the equation 


дх.` 


2 


== wet 


E Е а е 
| x: 8x; 9 М, M, tk + 5 М; а y (2;, ш) = 0. (14) 


which leads to the spectrum 


:— fat — +1) 


where it should be assumed that for the space-like vector и, а = 0 (this case will 
always beimplied below). For В > 0 and ғ > 0 the spectrum (15) has a point of con- 
densation for то = 0. If, for example, f > 0, but = < 0, then there will be solutions 
with my? < 0, corresponding to functions of the type V = е+= "4", Such solutions 
should, of course, be discarded. In [7] however, a more fundamental conclusion of the 
. mecessity of neglecting the equations leading to imaginary masses was given (we shall 
return to this in § 3). Thus we abandoned further investigation of equations of type 
(14), since it turned out that they always had either a decreasing spectrum or allowed 


E,? == тұ? == (15) 


solutions with imaginary masses?. 
The relativistic equations for scalar or bispinor functions Ч (x;, и) which were 
subsequently examined in a number of other papers [8—14] also possess the charac- 


2 Included in this group are the equations considered in [7], namely 


з 
(в — +byi yr М 2 V (xi, uj) = 0, 


' Эх; 
9? T b M. 77 29 ауан 
‚ 9%, Oxy ааны "Құла 7 


where, іп the first case, YZ із а bispinor, and, in the second, а 4-vector. 
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teristic peculiarities noted above or have solutions corresponding to zero rest mass. 
As we have already pointed out in [15], the well known equations introduced by Yukawa 
[8—9] are entirely of type (2). Indeed, the following equations were considered in [8] 


д 


| д. рош =, {шщ = ГУР —0, ur —0 (16) 


| x; 8x; 
It is quite obvious, from what was said above, that the system (16) corresponds to 
a particle with an indeterminate spin and not to a group of particles with spin zero, 
as, at first, assumed in [8]. In [9], the following equition was discussed 


92 u? 32 1 n 
xo 0x Ст ами (ани. = 0 17 
С дх; 2 | ди; ди; T us ши (x и;) ( ) 


9x; 


where in [9], м was denoted by А and и, by r, 
Making the substitution V^ (x;, uj) = v (ж) p (ш) we are lead to the solution 


9 (ш) = Hm (ulu) * Низ (изм) * Низ (Unal 2) * Hno (иуи) X 


2 2 QA 
X exp | au аан (18) 
2u 
ETY | 
тоё = -z (пү- ng + пу — т + 1), Пол„з›з = 0,1,2,.... 
u? 


where и, = iu, and Н, 


is a Hermite polynomial. 
Spectrum (18) is infinitely degenerate and contains values of zero for the rest 
mass. In order to remove the degeneracy, in [9], an equation containing mixed terms? 


was considered, namely 
| 92 A 22 1 2 
[95 dx, Уор ч ey sit ae a 


озы ا‎ 2 
+e |- (5) + 2 (u ze) || v=o 09) 


Solving this equation in the rest system, we obtain the mass spectrum 


ee 
xi + -a (п F na +t tg = ng + D? 
т? сала (20) 


1-26 (a 44 


According to the sign of e, as in the case of spectrum (15), this spectrum either 
contains a decreasing branch (ту — 0) or leads to the appearance of imaginary 
masses. i 
C cc A M T capu EC a ee c. us 


3 [n equation (19) of [9], е = В? and x? = 0. We introduce the value x? Æ біп order to get rid of 
a solution of the type V (x, w) = Ф (ш). 
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In [10], the following system of equations was used 


д д 
ЕЕ l S Hg m 
б x; 2); ди; == x) P (Xi ш) 0 


22 
Эх; 2x; 


и; 


вто, (в 


9х; 


Е) ф = (21) 


which, however, is not a consistent system. 

There is no need to dwell in detail on this point or on other critical observations 
made by us earlier [16] in relation to Ref. [10] sin e their basic features were subse- 
quently discussed by Rayski himself in [11]. 

In [11], the following system of equations is examined 


әз Е 3 ES әз 
-A сұғы =4¢| - ди; ди; tuu | э, n ди; ди» alg 


se Із ШЕ оғ Сал іе E а-а, 212) у 
58 ди; A ама ды дщ; PR Эх; 9x; HO ғ ( ) 


This system is a consistent one and leads to the spectrum 
1= 0,1.2,... (23) 


Неге we again meet with imaginary masses аз well as with differential supple- 
mentary conditions. 
The following equation was suggested in [13] 


a 9 
9? т Ox; дхь 
2. а жі зі Born ee | 4: 
Әл; Әл; ТЕ; P anc one 
9х; Эх; 


where M,, is the operator (3). 

An analogous equation containing the operator 93/2, 9x; was also suggested in 
[12, 14]. In such cases, we are no longer dealing with second order differential equa- 
tions with respect to x;, to which we are restricting ourselves. We shall note only that 
if both parts of equation (24) are multiplied by 9*/9х; 9x;, then we obtain a 4th 
order equation, two solutions of which correspond to zero rest mass. That is why 


4 [n [10] as well as in the analogous case in [12], an attempt was made to examine the supplemen- 
tary conditions of the type и; (9/92;) ¥ == 0 not as a new equation for V", but as a condition for limiting 
the region of variability of the variable и, (for example, from и, == 0, it was concluded that (0/u,) Y^ = 0. 
It is not obvious to us if such an approach is, in general, permissible in relativistic theory, and in any 
case, there is no doubt that such a procedure extends outside the framework of the general scheme [4] 
and meets with difficulties when an interaction with the field is introduced. We shall not discuss such 
attempts in detail, arid we note only that with such an approach, it is obviously not allowed in any case 
to proceed from the system (21), which has a quite definite mathematical meaning and is not a consistent 
system. 


170 V. L. Ginzburg 


even a transition to a 4th order equation, which is still not a particularly drastic 
change (it is simply associated with the doubling of the number of states), is not 
of interest to us in the given case if we are not dealing with particles having a zero 
rest mass. As regards equations of type (24), which contain derivatives in the denomi- 
nator, it is not clear to us, in general, how to work with them. 

Finally, we will consider the system of equation suggested in [14] (in (25) У is 
a hispinor) 


2 


д 
P $c val == Sud tuu) res и; = 0 


E! («- 2 )у-0 (25) 


This system is a consistent one and leads to the spectrum 


My = + [x + a (т + п + ng --1)), (26) 


which does not possess any of the ,,defects“ of the imaginary or zero mass type. 

§ 3. All the equations introduced possess characteristic peculiarities which 
ме already mentioned in discussing Ref. [7]. These particuliarities are as follows: 

If we have one equation of type (2) in which there are no mixed terms (as regards 
variables x, and u,), then the spectrum will be infinitely degenerate (see (6) — (11), (17)). 
The case of one equation of type (2) is then related to the general scheme [4] based 
on equations of class (1). . 

If we have several equations of type (2) relating to one function W being differen- 
tial equations with: respect to x;, then the system is overdefined and extends outside 
the framework of the scheme [4]. At the same time, such a system can have non trivial 
solutions in which the degeneracy, generally speaking, disappears. The systems of 
equations (6) and (12), (22) and (25) are examples of this; however, in the case of 
system (16), the degeneracy is preserved. As in § 2, we will call systems of two or more 
equations associated with one function ¥, equations will supplementary conditions, 
even though the basic equation and the supplementary condition sometimes differ 
from each other very little. Equations with supplementary conditions can possess 
a quite "good" spectrum, i. e., for example, a spectrum increasing to infinity with 
finite degeneracy (see (25)— (26)) and having neither an imaginary mass nor a rest 
mass close to zero. However, equations with supplementary conditions possess a con- 
siderable shortcoming — the field cannot be introduced in them by pos usual simple 
method. The situation here is similar to the case of particles with spins 7 3 or higher [5,6] 
for which it is not permissible to introduce the: electromagnetic field into the equa- 
tion by changing oes — d/dx, —ie Жы as it makes the system inconsistent. For | 


TS dn Me UA hereafter, to supplementary e we shall have i in ийла differential supple- 
mentary conditions, and not ‘conditions of the type (шш; — г?) y — 0. 
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particles with spins 2 2¢- 3) and so on, this difficulty сап be circumvented by 
introducing new functions and a variational formulation of the problem [5, 6, 1]. It 
is highly probable that an analogous approach is also possible for equations of type 
(2) with supplementary conditions of the same type. However, this problem has not 
yet been solved. At the same time, even if a solution to this problem is found, equa- 
tions with supplementary conditions will be very complicated if interactions are pre- 
sent. Therefore, the possibility of working without the supplementary conditions is 
especially attractive. 


The degeneracy can be removed not only by imposing supplementary conditions, 
but also by making a transition to equations of type (2) with mixed terms. In this 
case, if imaginary masses do not appear, then in all the known examples (see (14) — (15) 
and (19) — (20)), there appears a branch of the mass spectrum which tends to zero. 
This result has a very general character, as it is true [4], generally speaking, for the 
case of any equation of type (1), which also includes equations of type (2) with mixed 
terms but without supplementary conditions. It is thus clear that a non-degenerate 
increasing mass spectrum with one equation of type (2) can be obtained only in the 
presence of solutions with imaginary masses®. Solutions with imaginary masses 
V ~ exp {+ [mo] t) should, of course, be discarded, but it is not yet clear if this means 
that the corresponding equation itself is unsuitable. In [7], such equations were never- 
theless rejected, but in [9] only the inadmissible solutions were discarded, while the 
equations were retained. - - 


Up to the present time, in all applications to elementary particles, the imaginary 
masses do not appear. On the other hand, the rejection of one or another solution is 
not an unusual procedure in quantum theory, and, in particular, the discarding of solu- 
tions with imaginary masses can be treated simply as a generalization of the require- 
ment for the finiteness of the solutions. Indeed, the solution with E? = т% < 0 
and momentum p = 0 under a transition to a moving system of reference corresponds 
to a solution which has an imaginary momentum and which does not remain finite 
in space. Therefore, by intfoducing an additional requirement, we come to discard 
only solutions with p = 0, which increase without limit as ¢ > + oo and which аге 
obtained by a Lorentz transformation from inadmissible solutions with imaginary 
momenta. By virtue of the above, the discarding of equations with imaginary masses 
in [7,4] does not seem well-founded. However, one cannot also have full confidence 
in using such Шш since the presence of the discarded imaginary solutions can 
“take revenge“ when we deal with the interaction of the particle with the field. In ad- 
dition to this feature, to which we shall still return, it should be noted that the presence 


> د د‎ ee 


4 We note that the statement about the appearance of a decreasing spectrum for equations of type (1) 
represents the rule, but by way of exception an increasing spectrum сап also be obtained (see, for example, 
(9) and (11)). This spectrum, however, is infinitely degenerate. It is probably not possible to obtain 
a non-degenerate increasing spectrum (without the presence of i imaginary iode à although this state- 

ment has. not been Jm proved. 
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of imaginary masses usually accompanies also the appearance of space-like solutions 


(solutions with space-like momenta p; [16]). 
Let us examine, for example, the equation 


9? E 
—— +3 ues 21 
к 9х; Tx ) р (х) ( ) 
which in the rest system has solutions with imaginary masses, since Е? == m? = — x?. 


We discard these solutions and other solutions with imaginary Ё and р as a result of 
applying the generalized principle of finiteness of solutions (see above), but this equa- 
tion will still have solutions with real Ё and р satisfying the relation 

Е? = — x? + p’. (28) 
Such solutions for which p> x? are everywhere finite and it is impossible to discard 
them on the grounds of the required finiteness. On the other hand, these space-like 
solutions correspond to velocities greater than the velocity of light (v = 9. E/9p = 
= р[\ p? — х? >1) and are obviously inadmissible because of physical considerations. 
The simultaneous appearance of imaginary masses and space-like solutions is the rule, 
and it occurs, perhaps, for all equations having a finite number of dimensions (this 
problem has not been investigated in detail). However, for equations with an infinite 
number of dimensions, cases exist where there are solutions corresponding to imaginary 
masses, but no space-like solutions. For example, equation (19) for E = 0, p, = р» = 
= 0 and p, >20, has the solution 


2 
DIEN VA pU e oco (29) 


If e > 0, equation (19) has solutions with imaginary masses (see (20)), but in this case 
p^, <0 everywhere and therefore there are no space-like solutions for which p? 
would be greater than zero and E equal to zero. On the orther hand, for many other 
equations, space-like solutions appear. This occurs, for example, if in (19) the third 


term is equal to 
и? 9? 1 
апд зо оп. | m E дщ Тя М u) 

As shown in [17] space-like solutions also appear for the system of 
equations (22) considered by Rayski [11]; in this case there are also solutions with 
zero mass. All similar equations containing space-like solutions apparently should 
be considered inadmissible. I é 

As regards equations with imaginary masses, but without space-like solutions, 
there is, as already pointed out, no basis for discarding them as long as we are dealing 
with free equations. Thus some equations with imaginary masses and also individual 
equations with supplementary conditions obviously lead to increasing mass spectra 
with non-infinite degeneracy. In addition, one can apparently propose as many such 
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equations аз one wants by modifying the corresponding form of function Ф in (2). 
For these reasons, we somewhat sceptically regard the attempts made [11] to connect 
one or another equation of type (2) with the mass spectrum of known particles; if 
we are not to rely simply on a lucky guess for choosing some equation out of all possible 
equations, then we need a new idea which so far has not been enunciated?. 

Despite the fact that it is difficult nowadays to connect the relativistic wave 
equations with mass spectra with well-founded hopes for progress in the sphere of 
the physics of elementary particles, further investigation of these equations would be 
of interest, if even from methodical considerations. 

Thus the centre of attention undoubtedly should be shifted to the examination 
of equations with interaction and investigation of the solutions of those equations, even 
though in the simpler cases (constant external magnetic field, onedimensional electric 
field, first-order radiation processes, and so on). In speaking of interaction, we have 
in mind here the usual local interaction, the study of which, one would think, should 
precede the introduction of the non-local interaction, since the examination of the latter 
is accompanied by additional difficulties. 

In our opinion the above remarks apply in even greater measure to the case in 
which one considers infinitely degenerate equations in the hope that the interaction 
with the field will remove the degeneracy [9]. In such conditions, perturbation theory 
is not applicable (as it cannot, generally speaking, lead to the appearance of new values 
of mass) and the interaction problem should be solved by more accurate methods. 
But the application of such methods, as it is known, meets with great difficulty even 
in the case of local interaction with the simplest fields and the transfer of these methods 
to the equations interesting us is by itself an independent, very difficult problem even 
for local interaction. Finally, we note that the presence of strong interactions makes it 
possible to use the free equations for finding the mass spectrum only under the assum- 
ption that the change of mass due to the interaction vanishes or, at any rate, is very 
small. Such an assumption, however, is very unlikely if we are dealing with a theory 
without divergencies where there are no grounds for renormalization, or simply leaving 
out of account all the field contributions to the self-energy of the particle. Thus, 


7 [f new tensor arguments of the type S;, are introduced into function У, then, generally speaking, 
a continuous mass spectrum is obtained. By considering V/ to be a function of some scalar, depending 
2 Ж 
on the form of function Ф in the equation — —x? + ø (^ x| |, (xj, u) = 0, апу 
9х; Ox; eu 
spectrum can be obtained. The introduction of a scalar argument is particularly simple and does not 
cause any difficulties, but on the other hand, it can hardly be of any interest in connection with - 
physical applications. 

8 We are now making no mention of the fact that in recent times there has been no convincing 
basis for taking the spin of mesons and hyperons to be larger than 0 and 1/,. If the spins of all particles 
are equal to zero or 1/,, the mass spectra probably cannot be obtained by considering variables of the 
type uj, which although not spin variables in the literal sense of the word are clearly connected with the 


spin of the particle. 
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the possibility of comparing the spectra of the free equations with the mass spectra of 
real particles becomes still more illusory. 

We shall leave, however, the discussion at this point, since the principal purpose of 
the present paper is only to cast some light on the present state of the problem concerning 
the properties and character of free relativistic wave equations with mass spectra. 


KPATKOE СОДЕРЖАНИЕ 


В. Л. Гинзбург, О релятивистских волновых уравнениях со спектром масс. 


В последние годы в ряде работ обсуждаются релятивистские волновые уравнения со 
спектром масс, содержащие новые ‚непрерывные независимые переменные. При этом, 
однако, не обращается обычно внимания на ряд затруднений и особенностей, связанных 
с исползованием таких уравнений. 

Релятивистские волновые уравнения, зависящие от новых переменных, в большинстве 
случаев от некоторого пространственно-подобного четырехмерного вектора u;, подробно 
рассматривались в работе [7], а затем и в ряде других статей [8—17]. Таковы, в частности, 
уравнения (2), (6), (14), (16), (17), (19), (22), (24) и (25), приведенные в тексте. Все эти урав- 
нения имеют следующие особенности: 

Если имеется одно уравнение типа (2), причем в нем нет членов смешанных относително 
переменных =; и Uj то спектр масс бесконечно вырожден (см. (6) — (11) и (17). Случай 
одного уравнения типа (2) относится при этом к общей схеме [4] базирующейся на уравне- 
ниях класса (I). , 

Если имеется несколько дифференциальных уравнений типа (2) относительно одной 
функции V (a; uj); TO система выходит за рамки схемы [4] и, вообще говоря, переопре- 
делена. Однако, подобные системы могут иметь нетривиальные решения, причем вырож- 
дение обычно снимается (см. системы (6) и (12), (22) и (25)). Основной недостаток этих 
уравнений, которые мы называем уравнениями с дополнительными условиями, связан 
с трудностью ввести в них взаимодействие с другими полями. 

Если имеется одно уравнение типа (2), содержащее смешанные по хи, члены, вы- 
рождение спектра также снимается. При этом, если не появляются мнимые массы, то во 
всех известных случаях имеется стремяшаяся к нулю ветвь значений массы покоя (см. 
(14 — (15) и (19) — (20). Поэтому получение невырожденного растущего спектра масс, на 
базе одного уравнения типа (2), связано с наличием решений с мнимыми массами. Эти 
решения должны, конечно, быть отброшены, что еще быть может не делает непригодным 
само соответствующее уравнение. Однако, в целом ряде случаев появление мнимых решений 
оказывается связанным с появлением также и пространственно-подобных решений (ре- 
шений с пространственно-подобным импульсом pj. Эти решения физически недопустимы, 
а отбросить их с помощью принципа конечности решений нельзя. Поэтому сооветствующие 
уравнения, к числу которых принадлежат уравнения (22) из [11], повидимому, должны 
заведомо считаться недопустимыми. 

Таким образом только некоторые уравнения с дополнительными условиями и уравнения 
с мнимыми массами, но без пространственноподобных решений, не вызывают возражений 
общего характера, пока речь идет только о свободных уравнениях. Не исключено, однако, 
что при наличии взаимодействия попытка использования и таких „хороших“ уравнений 
столкнется с новыми серьезными затруднениями. Исследование свойств и решений реля- 
тивистских волновых ‘уравнений со спектром масс при наличии взаимодействия с электро- 
магнитным и другими полями представляется нам имеющим несомненный методический 
интерес. Что же касается использования подобных уравнений в области физики элемен-. 
тарных частиц, то в настоящее время трудно надеяться на успех в этом отношении. 


[жип 
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SOME REMARKS ON THE THEORY OF THE 
ELECTRON — PHOTON CASCADE 


By Jan LoPuszaANsKi 


Institute of Theoretical Physics of the Polish Academy of Sciences, Wroclaw 
(Received February 16, 1956) 


The well known integro-differential equations for the first moments of the 
distribution for the electron — photon cascade are derived from a integral equation 
analogous to the Smoluchowski equation from the theory of stochastic processes. Some 
asymptotic properties of the distribution function are given. 


1. Introduction 


This note contains some preliminary considerations, which are useful for the two 
papers which follow (Stachowiak 1956, Konwent and Lopuszariski 1956). We give 
the derivation of well known integro-differential equations for the first moments of 
the electron — photon cascade in approximation A. In this derivation, we start from 
an integral equation, analoguous to the Smoluchowski equation! in the theory of sto- 
chastic processes. In addition we present some remarks on the behaviour of the distri- 
bution function. 


2. Notation 
We denote by 
13 
— sik) l 
ne (e, x) (1) 
the average density of particles of type К in the energy range (e, е + de) E, at the depth 
x due to a primary particle of type i and energy Е’. The supperscrips i; ko РА 
where we designate the electron by 1 and photon — by 2. 

We have further 
' UE p í 
8% № (e, x) = — — — (e, x) (2) 

where 5% ® (e, x) denotes the average number 2 particles of type k of energy larger 
than е E, at the depthyx due to primary particle of type i of energy Ey. 


1 This equation is also called the Chapman — Kolmogorow equation. 
(177) 
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3. Derivation of the equations 


The starting point of our considerations is the integral equation 


1 
2 
қ E dé’ - 
g(t, k) (e, x) = | У 500 (г х) 50 5 (3 ‚® =) TT (3) 
6 1-і 


forix cor XU. 


We write the following obvious relations 


s% P(e, 0) = by 8 (1 — 8) (4a) 
— 6(1 = е) a(D + w (1 — 8) for J = 1 
b их) for] = lk = 2 
(% 580 № (в, 2 n (=) ог (Ab) 
9х x-0 2w (=) for L= k = 1 


— 0 (1.— е) о forl=k=2 


where 4? (ғ) апа 10% (=) denote the bremsstrahlung and pair production ditferential 


cross sections respectively and 
1 1 
a — f w (6) ав, | a = f w (е) de. 
б 0 


We differentiate both sides of equation (3) with respect to x and put ж = x. 
Taking into account relations (4), we obtain 


1 


i 1) 
95%5 (e, x) = — al sD (e, x) + [= D (e', x) x 


9х 
1 
x wd n M | = m9 fe 2) (г, x) w (=) = | (5a) 
б 
1 
956. 2 (e, 
a x) = — aO st D (e, x) + | 50 р (e’, x) wt (=) = (5b) 


0 


This is already the first group of well known equtions, which are discussed for the case 
. of boundary conditions (4a) e.g. by Landau and Rumer (1938), Bhabha (1950), Bhabha 
and Chakrabarty (1943, 1948). 

To obtain the second group of equations, we differentiate (3) with respect to x’. 
Further, we make use of equation (3), differentiated with respect to x, and finally 
рих c0. ; 
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With help of (4b), we have 


9 sb B (siu) й 
2x ж 


1 
чу ] ә (e) le k) ( ы a :) +. 502, 9 (5 ; Jl = (6a) 
0 


1 


— a (D sa k) (ғ, x) 4 


апа 


9 s( 0 (e, х) 


5 = — а 50 9 (e, х) + 2 | w?) (є') 506 ® | 
х 


0 


E dé’ 
Улу 15 x) P (6b) 


€ 


These equations form the second group of integrodifferential equations. 
Integrating (6) into = т the interval (у < 1,1) we arrive at the well known equations 


1 
950 ® (m, x) + al) SQ. № (п, x) = fr (е) fso М ( : а) I 
c {2 
0 


Эх 


+ 50,5 (2 :)| de’ (7a) 


€ 


and 
1 


a® SG. ® (m, x) = 2 Г их?) (e) 50 ® (2. de’ — (1b) 


0 


25%, ® (1, 
— i 
x 


considered e.g. by Janossy and Messel (1950). 
In the first paper following this note extensive use is made of (7). The second 
paper makes extensive use of (5). 


4. Behaviour of the distribution function 


We shall now show that 5% P (e, x) is of the same order of magnitude as the proba- 
bility of finding at least one particle of type k at depth x and having an energy larger 
than e E, along with an arbitrary number of particles of type (3—X) due to primary 
particle of type i and energy Ey. We denote this probability by PC ? (e, х). To show 
this, we note that 


B SG 9 < pa» (e, x) < St Ne, x) (8) 
€ 
Urbanik (1955) has shown that 
| Sa k) (e, x) = е-а(2)х ye Ё) (є, x) onl (9а) 
where for each В > 0 , | 
lim е % }7® 9 (в, x) = 0 (9b) 


X00 


180 J. Łopuszański 


and 
lim WE 9 (e, х) = оо (9с) 


Xx + со 


Let us write P^ ? in the form 


p ® — ега» VO ® (e, x) (10a) 
Then, from (8) it follows that 
lim e * V€ ® < lim e Wo” — 0 (10b) 
x ¬+ co x со 
and further that 
J-i 
lim V@ № > H lim Wt ® = со (10c) 
± +00 € х-со 


Equations (10) represent the generalization of Urbanik’s result for the distribution 
function. 


5. Final remarks 


From (6) and (7) it is obvious that for large depth of the absorber all asymptotic 
properties deduced from (7) for 5% 9 are also valid, mutatis mutandis, for 5% 9, 
Further, from (8), it follows that the same is also true for Р 9, 

Consequently the results obtained in the first of the two following papers (which 
deal with equations (7)) are also valid mutatis mutandis for 5% ® as well as for P 9, 
In the case of the second paper, which deals with equations (5), a similar situation 
occurs, but further proof is needed. 


KPATKOE СОДЕРЖАНИЕ 


I. Лопушанский, Несколько замечаний о теории электроно-фотоннго каскада 


Из интегрального уравнения, аналогичного с уравнением Смолюховского в теории сто- 
хастических процессов, выведены известные дифференциально-интегральные уравнения 
для первых статистических моментов для электроно-фотонного каскада. 


Кроме того представлены некотороие асимптоматичные свойства функции распределе- 
пия для болыпих глубин абсорбера. 
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SOME PROPERTIES OF DISTRIBUTION FUNCTIONS 
OF ELECTRON — PHOTON CASCADES 
AT LARGE ABSORBER DEPTH 


By Henryk STACHOWIAK 
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It is shown that at large depth of absorber, the number of particles found in elec- 
tron-photon cascades depends only on the spectrum of the primary photons, and that 
at large depths, the number of particles decreases rapidly as the energy of the particles 


increases. 


1. Introduction 


The object of this paper is to investigate certain asymptotic relations between 
first moments and quantities connected with them which characterize the probability 
distribution for an electron-photon cascade at large absorber depth for threshold 
energies higher than the critical energy, since the work in this paper is done in appro- 
ximation A. We have not made any assumptions as to the functional form of the cross 
sections except that we require the cross sections to be bounded functions which do 
not vanish within a finite interval in the range of their arguments. We further require 
that the functions will be homogeneous functions of the energy before and after colli- 
sion (this condition is satisfied in the ultrarelativistic approximation with complete 
screening) and that the cross section for pair production will be smaller than the cross 
section for ,,Bremsstrahlung“. 

A similar treatment for nucleon cascades in a homogeneous nuclear matter was 
made by Urbanik (1955a) and Łopuszański (1955). The computational D used 
in this paper are different from those used by these authors. 


2. Notation 


The electron-photon cascade is described by the function P (nj, ng | е, x) where 
= 1 signifies that the primary particle was an electron andi = 2а photon; PO (пу, л, | 
е,х) is the probability of finding n, electrons and ng photons of energy > gE, in a cas- 


1 We assume, that the primary energy spectrum is monochromatic i.e. that all primary 


particles have the same energy. 
(181) 
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оь OSA лайым мыл ТЫСЫ ЕЕ: а 1. у... зе ыш 
cade at a depth х in the absorber, if the primary particle was of the type т, where Ey is 
the energy of the primary particle. 

If we are interested only in the number of photons. or only in the number of 
electrons, we introduce the function P® ® (n | e, x) which represents the probability of 
finding n particles of type A. 

We obtain the relation 


oo 
PEM (nile x) = № P (n, , n, | в, x) (1) 
п3-ь-0 
The avarage number of particles is given by the formula 


SOP бы д = 5 n PO 9 (n | e, x) (2) 


n=0 


We introduce the generating function 


СӘ (e, u, Ua, x) = р. > РФ (y, n, |=, x) иие (3) 


=0 m=0 


Function GP (e, и, x) is defined as follows: 


Сели) = СӘ (e, и, 1, x) POUT | e, x) u” (4a) 
n=0 
GÈ (e, u, а) = GP (e, 1, u, a) = У PÈP (n | e 2) и" (4b) 
| n=0 
We note that 
Ge 9 (e, 0, x) = PÈR (0, е, х) (5) 
апа therefore 
PED (> 0| г, х) = 1L— СО (0, х). (6) 


Р®® (50 | e, x) represents the probability of finding at least one particle of type k with 
an energy > е Еу at a depth of x. For simplicity we will write Р®® (>0 | ғ, x) as 
p bg x). 

The Jánossy equation (Jánossy 1950) holds for the ВОЕН function in appro- 
ximation A. Thus. 


IG® s Uls Й ; 
209 e و‎ + а; GO (e, uy, Ug, x) = 


1 
- [eo Ба a Uy, Ug, X x) GG-) E uae tas x} ш (ғ?) de’ (7) 


0. 


where 29 (e^) is the differential cross section for a particle of type i. 


1 
ge © (e) de’ | 
а ІС (е) de (8) 
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3. Behaviour of the Mean Number of Particles for Large Absorber Depths 


We will show that the following relations for the mean number of particles held to 
a first approximation: 
Ех 
ag С 
х-оо SG ) (е, x) 


9059 (e. x) 


= 0 юге >e (9) 


їп = 10 

ae 50.9 (в, x) ou, 
(2,9 (ғ” 

т eer 2212) DOS T (11) 


хо 500 (ex) | 


H 


Е :) | 
lie, 12 
lim (0% (ғ) de’ ee ау — а; з, 


= 50 (s, х) 


Proof: 
As the starting point of our proof we will take the well-known equation for the 


first moments (Jánossy and Messel 1950) 
1 


25386, а, k) z ЕТЕ Кезі 
hasc 2 ` a; S (2, x) = 5 МЕСТЕ Ned + 


0 


4-58-20 (5 - 3! wf? (e) de’ (13) 


К. Urbanik (Urbanik 1955b) showed that the function 5% ® (e, x) satisfies the following 


relations 


lim e*m * S9 Ю(є, x) == oo а, = min (оу, a5) 
x со 

lim es^! e°m* St ® (e, x) =0 8 >0 

x- oo . 


As a result, е%* 5% 9 (e, x). increases more slowly than an exponential function. 
We thus can write 
SE (e, x)= етн" Д/% P (e, x): (14) 
where W ® (е, x) increases at a less-than-exponential rate. Е 
Inserting equation (14) into the left-hand side of equation (13), we obtain 
OWE № (е, х) _ 
Ox EE 


1 
- f ш (е) de Е »( : а) + 38-59 кү (15) 
is | € 


(a; — Am) ет: WE 9 (e, x) + ет 
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From collision theory, we know that a, > аҙ. Thus а, = аз. Dividing both sides by 
S@® (e, x), we obtain 
1 IWG 0 (е, х) | 


1 
| ба, ю s Se :) а-ы» (2. б 4| 116) 
05 РОТ E RE an ai T; 
«9 (e) de 56 ® (e, x) 5% (e, x) 
0 


We note that since W“ ? (e, x) does not increase as rapidly as an exponential fun- 


ction, 


1 OW ® (e, x) 


e гету игта о 
Therefore, as a result of (16), we obtain 
Sa, ® s E 7 SG- i, ® Es x} 
i i) (= d © у N A r ELSE AL) See 
lim, Барма a S. 9 (e, x) о (е, x) а— (17) 
0 | 
Setting i = 2, and noting that w® (e’) = wu? (1—'), we have 
sa. e( 5, i) 3 
537 ая аа ЖТ 


Since 56 ® (e, x) is a non-negative function and w (e) > 0 within the limi ie of m 
integration. 2e | 
: DAE A; 


In order to ерда 
Уш: 


4- 


м 


DA ا‎ Аза 
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For т = l, equation (17) has the following form: 


1 
у є € 
Sa, »(; 6 x) SQ. k) (s 5 x) 
lim wD) (e) de’ E 5 


—— On — 
А 7 $0.8 (s, x) Sd. ® (e, x) NE қ. 
0 


(20) 


If we insert into (20) the two relations obtained previously, we obtain a result 
contrary to the assumption made in cascade theory that w® (&') is bounded. Thus the 
2, 4505) (exa aa ge қ ; 
assumption that lim “°” is finite turns out to be incorrect. This ratio tends 

wees OU 2 esu) 
to infinity. 
59 ® (e, x) 


Therefore (10) is correct: ISO т 2:0 


лб. x) Sq. k (5 | ») 
eee ё ; E 
Multiplying tm, 500 (ea) by m. ume AES 
senf, x), 
ё 


we obtain an expression analogous to the first term in the integrand of (20). This 
means that the first term in the integrand is small in comparison with the second, and 
therefore can be neglected. We thus obtain (12): 


1 
1 € 
50, Ю 2 р x) 
lim GO (e y nmen ee 
= (е) 90 9 (e, x) re 
0 


Also, 


sa, ТЕ Ж >) 
lim wD (e') de’ ағы UE ев 


× ¬+ со Sa k) (=, х) 
0 


From this last equation, we obtain 


50,9 (ғ, x) 
im => esce e) fore е 
GUSTO (ea) 


which is the first part of (9). Proceeding in an analogous manner, we multiply (12) 


(1, 0) 
by Jm pe and obtain the second part of (9): 
б 5,0) (ғ, х) / 
um 50,9 (s, x) 0 for е LE 
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It now remains for us to prove Eq. (11). To do this, we will write Eq. (12) in the 


j 1 
50, (6 1 =) е 50, (5 " г) 
sto) p] ^w (9) ded — e 


lim wD (e) de’ “SGP ТЕЙТ" SGP (e, x) 


form 


„ 


(up c 02- 
x +00 


where (1-е”) is a finite number. Because of (9), the first integral with the argument 


Ox 
— > = can be ignored in comparison with the second. 


Thus 
B 4 
50, Ё) А us) E 
e. 
lim wD (e) de © ( 
x +00 (e) SLA (e, x) \ 
0 d 
as a result of which | : А b3 "A 
SQ. 5 (e. x | AD 4 1 
Г for & > Е 


HW LÀ 


4. Extension of the Obtained Results on Р®® ferc 


Because of the relation (see the preceding paper, Łopuszański 1956) 


LIA 2. M. Em 
В SG 9 (e, x) ei (e, x) е х). = | | bee 


the ratio of 5%? (e, x) to P69(z, x) ва certain. finite 5 ТЕ ка than unity. 
Therefore, from Eqs. i (10), ШЫ? and (12), we obtain | - 
= 


1 


E ACO (eran) == һө (2) de’ 
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рк pe 
| 1—ғ'° 
0 
+ pa-i k) (5 ; x) _ Ра,» ( € ге ) . РЗ-БЮ (5. 3] (26) 
= = 6 


Since from (14) and (21), P®® (e, x) decreases exponentially with the depth, we 
can, in the limit, neglect the last term in the integrand of (26) as being of small magni- 
tude to a higher order. For large x, the equation 


9 PG. 9 (e, x) 
Ox 


+a;P@ (e, x) = 


1 


= [ (€) dé’ | Po k) (25 J =) == pou a (5. 3] (27) 
0 


is satisfied This equation is analogous to equation (13). 

Because of (21), Р® P(e, x) can be written in the form PC 9 (e, x) = e-*m * (е, x), 
and therefore, the same relations are satisfied for Рё ® (e, x) as for S6? (e, x). Thus (25) 
takes the form 


1 


po» (2 ) 
7209 8 
lim s (ғ? de’ уа АГ ины дшге @ — 0%. (28) 


PGP (e, x) 


x> со 


5. Some Properties of the Conditional Moments for Large Depths of 
Absorbers 


If Po» (п | £, x) represents the probability of finding n particles, if there is at least 
one particle. Then the relation 


P9 (п | г, x) = PO” (n | e, x) РФ» (> О | e, x) for n 0 (29) 
will be valid. 


We define the conditional moment (the mean number of particles if there is at 
со 


least one particle) to be “SGP (сул) = пРФ%9(, Г) (30) 
= 


Because of (2), (6), and (29) 
S6. (e, x) 


SG k) (e, x) == P9 (s. x) [5] (е, x) ae 


(31) 


25% № 
Evaluating Да, 


2 on the basis of Eqs. (31), (13), and (26), we obtain the follo- 


wing equation for SC? (е, x) 
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en | d >») 
95% (е, х) == we) (2) de’ ul 21 (5% o(; : 4 x) MEA (e; J 3 
t 


Dx ро Ё) (ғ, ж 


0 
1 


pG- i» 19 à x) 
= : € TU 
aero IET 


0 


1 


€ с € 
› k) (3 — ty k) — 
Ж | E = ' =) k (: J (32) 
— 509 (в, x) wh) (a) det ЕНЕ ые 


P& ® (e, x) 


0 


хы. 50509 
We assume that 5% 9 (e, x) remains finite, then lim 
х» со x 


— 0. [n the limit, we 


can ignore the last integral in (32). We then obtain 
1 


| ё i Е | 
А 25679 T 
ae ке Р®% (s, x) E (=; : >) е J 


x ¬+ оо 


1 


рз-ью (> : >) 
кш) و چچچ مه(‎ [rt (oss) Зена =0 69 
0 


Because of (22) and (23), we obtain an identity be letting i = 2. If we let i = 1, the 
first integral vanishes by (22), and there remains 


1 


pan (=>) (34) 
3 ; р € = € = < 
шп и (=’). de =P Paes E k) (5 қ ) - 80,9 (е, J =0 


х» со 
0 


Using the mean value theorem, we obtain, because of (28) 


(a, — dy) lim Е ® ( А s) — 50,8 (e, » | =0 (35) 


End e” (x) 


Because of (24), we can use the small interval (e’’, 1) for the limits of the integral 
of (34) instead of the limits from 0 to 1; ғ” (x) will then lie near 1. 
Since 5%» (e, x) is continuous and bounded we can pass to the limits 
lim [5 (e, x) — 509 (в, х)] = 0 (36) 


х- oo 
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Since $®® (е, х) is always > 1, Eq. (36) signifies that in the limit 
SO (e, x) = 509 (e, x) (37) 
In the paper which follows (Konwent and Lopuszaüski 1956), it is shown that 


lim S69 (e, л) = 1 


х оо 


6. Physical Conclusions 


1) The spectrum of particles at large depths is very steep, i. e., if we somewhat 
raise the threshold energy, considerably fewer particles will be observed at large depths 
of the absorber. 

2) If, at a large absorber depth, a particle is encountered, and if the primary parti- 
cles are both photons and electrons, then the encountered particle is due to a photon; 
the electron-photon cascade does not display any ergadic properties. Only the primary 
photon spectrum influences the number of high energy particles counted at large depths 
(see footnote on page 181). We thus see that the ,,privileged“ photons exist not only 
in the approximation taking into account ionization losses and the penetrating power 
of low-energy photons, but already exist in the ultra-relativistic approximation A. 

The difficulties of experimental observation of this effect results from the fact 
that the ratio of the probability of encountering a particle due to a photon to the pro- 
bability of encountering a particle due to an electron increases more slowly than an 
exponential function. while the probability of encountering a particle decreases 
exponentially. In order to check whether or not the effect has been observed, numerical 
calculations would have to be made. Thus the problem is to decide what is to be consi- 
dered a large depth for a given threshold energy. At a depth where the above-mentioned 
ratio of probabilities is of such a size that the primary electrons can be neglected, the 
density of counted particles will be too small if the ratio increases too slowly. However 
an increase of this ratio with an increase in depth would confirm the results obtained 
in this paper. 

The results we have obtained are weakened by the faci that the spectrum of 
particles at large depth is very steep, and if we somewhat raise the threshold energy 
of the particles derived from protons, and keep the energy of particles derived from 
electrons constant, the ratio of probabilities will decrease to zero with an increase in 


depth. 
ра.» (> ШЕ J 
» ^ о (1) , CO Up Gea iow. 
3) The functions Ug w™ (e ) Pa. (> 0| в, x) 


SQ. ® (5. :) 
є А 
і (D (gy — ———À-—— ^ are of the Dirac д type. 
and E «t? (e) Sd. 9 (e, x) УР 
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4) If, at large depths, we take into account only such cascades which have not 
as yet died out, then the mean number of observed particles will be the same regardless 
of wherther the primary particle was a photon or an electron. 

Note: The term, ,,threshold energy“ as used here should not be confused with 
the threshold of equipment since, in general, this is very low, while the theory employed 
here is valid only for high-energy particles. 

I would like to thank Doctor J. Lopuszariski for suggesting the topic, and for 
his valuable aid. I would also like to thank Professor J. Rzewuski, Professor M. Mie- 
sowicz and his co-workers, and Professor A. Zawadzki for their helpful comments. 


KPATKOE СОДЕРЖАНИЕ 


Г. Стаховяк, О некоторых свойствах функции pacnpedenenua электронно-фотонноео 
каскада на большой глубине абсорбента. 


Доказано, что для электронно-фотонной составляющей космического излучения число 
частиц с энергией>><Е, встречающихся под толстым слоем абсорбента зависит единственно 


регистрируемых частиц всё скорее уменьшается с возрастанием наименьшей энергии, KOTO- 


| от первичного спектра фотонов!, а также, что c возрастанием глубины абсорбента число | 


UU ма ал, АА 


| рой должны обладать эти частицы, чтобы они были регистрированы (наименьшая энергия - 


остается болыпе критической). 
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It is shown that 
т SG, 2) (e', x) 
in = 2 ee for e > г 
нь 56 D (s, x) 


апа 
SG. 2) (e, x) 
къо SH D(8, 2) — 

where 5% № (e, x) denote the first moments for the electron-photon cascade, and further 
that the probability distribution function Pf 9 (n, т | ғ, x) has the form 

PG(0,0| є, x) ~ 1 

P@ (0,1 | & x) ~ 50 2 (e, x) 

PG (n, m | в, x) = о (P (0,1 | в, x)) 
for (n, т) Æ (0,1) and n + m > 0. 


ТГ ntroduction 


The present paper is a continuation of the two preceding papers (Lopuszaniski 
1956, Stachowiak 1956, hereafter referred to as papers I and II гезр.). We make use 
of the result obtained there. We also employ the Urbanik’s conclusion (Urbanik 1955), 
that the first moments of the distribution may be written in form 


SEP (e, x) = eO WE (e, x) -- (1) 


(for the notation used іп this paper see paper I and 1), where Й7® (e, x) increases 
monotonically with increasing x and tends to infinity with increasing x not as rapidly 
as an exponential function. Further, we make use of an integral equation for the 


(191) 
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generating function obtained from the formal solution of the С equations of Janossy 
(1950) and given by Łopuszański (1955a). If we denote the generating function by 
С (€, цу, иә, x), this integral equation reads 


Г® (e, uy, ug, x) = (1— ыу) S€ D (e, x) + (1— и) SÈD (s, x) + 


x 1 1 
eae є' € € 
ON ET, d —— 561 (p, x — &) س‎ (5) Io Ет 7 гә 7 s Uy > =) = 
6 n| m 5 (x ) " E 1:79 n— 1:% 
0 0 0 


D 
€ 


Ж 1 1 
n Rg , j х ғ” E E 
E: : а | dn | 2 (n, х—ё) иФ A Га if С. Us, =] го = ЯТ ЕШ Ж s) 
б 0 (2) 


where 
1 (e, иу, из, €) = 1— СФ(г, шь us, 2) 


We consider the asymptotic behaviour of the distribution function and first mo- 
ments of the electron-photon cascade for large depth of absorber and positive threshold 
energy in approximation A. 

From the mathematical point of view, our treatment is not quite rigorous, e. g. 
we assume the integration or differentiation operators to commute with the limit 
operator. Some indications as regards this question are given in the earlier papers 
of one of the authors (Lopuszanski 1955), where analogous questions on the asymptotic 
behaviour of the nucleon cascade in homogeneous nuclear matter are considered. 

We assume that the sections for bremsstrahlung and pair production, w® (E, Eo) 


and w® (E, E) respectively, are homogeneous functions of the energies of the partic- 
les before and after collision, that is, 


~ 


wp vel 
MONT a oa eset oe 
( ? о) ES E, 
(in agreement with the well known expressions of Bethe and Heitler for the full screen- 
ing approximation). We also assume that the above-mentioned cross-sections have 


no singularities and finaly that in the interval O « E « E, no finite subintervals 

exist for which w ls is equal to zero. On the other hand no assumptions about the 
0 р 

shape of the cross section functions are made. 


We require, in accordance with field theory, only that 


1 1 
а® = f w® (e) de < a® = T w (e) de 
ò i 


. We emphasize, that from the inequality, given in paper I, for the probability 
ASE 


Pt 9) (e, x) < 56 0 (e, х) < B PG. № (e, x) | (3) 
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it follows, that the quantities Р% ® (e, x) are of the same order of magnitude as 5% 9(e, x) 
respectively (if P&P (e, x) tends to infinity or to zero, then SÈ ® (s, x) also tends to in- 
finity or to zero, respectively). Consequently all relations obtained for S? (&, x) 
also hold mutatis mutandis for РФ 9 (e, x). As seen from paper I, we can also write 


PEP (e, x) = ea x VOM (e, д), o 


where V9 (e, x) tends to infinity as x approaches infinity more slowly than an 
exponential function. 


Likewise we may easily generalize the result Urbanik to the functions 


ISEH (e, x) 


s6 № (e, x) = дє 
1. е. we may write 
5® (e, x) = е x wh” (e, x) (5a) 
where 
lim 206 9 (e, x) = oo (5 b) 
хә со 
а О (e, x) =0 B0 (5c) 
x + OG 


We note also that almost all relations for first moments, deduced in paper II, 
are also valid for 599 (e, x) because the equations for 5% ® (e, x) are very similar to 
that fulfilled by S? (e, x). In particular from 


569 (e', а) | 


аас 6 
ЕЛЕ РА (6) 
it follows that 
(7 k) (= 
lim — сар >e (7) 


SISTER D 


2. Asymptotic behaviour of some quantities characterizing the distribution 


The 5% ® (e, x) satisfy the following well known equations (see e. g. paper 1) 


(i, 1) f же; р E ағ” 
95% D (e, x) + aD 5% D (e, x) = fo D (e', x) wD (1 - A im S 
\ 0 


Ox 
1 E 
+2 [ و‎ ш® Se (8) 
ү aint 


€ 
1 = 
i, Р | de’ 
SIDE + a® 50 D (e, x) = ІІ 5% D (='. x) wu? (5) p (9) 
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Making use of (5a) we obtain from (8) and (9) respectively 
1 


Soter SOT REN Sem її 0% 2 (e', x) wO? | m | g ir 
Эх р е € 
0 
1 
а fosse aum (2) 
0 
| 
Iw (2,3) [cn лоо. [|4 
сш (е.а) шо (2) 4 ay 3 


Let us investigate the behaviour of these equations for x — оо. For this purpose 
we divide both sides of (10) by w® P (e", x), where ғ” < e. We get 


Iw; D (e, x) 


ожа А. 


tox ns Pate) (ery) Wie Denta) жүзү 
———Á шы (DE О ао CL) йе зел ее 
wes РЕ х) = (а а ) wE 1) (et. PS wh 1) Ee x) 1 =’ : 5' T 
1 ` 
wh D (e', х) de' | 
DID IM Е аа. да 
+ 2 J we» с x) их ^+ =’ у п) 
, 2 | TI 
From (5b) and (5с) it follows that w* (е, x) increase with x more slowly than an ех- | 
Ponential function. Consequently, we have = > 
90% 9 (в, x) | в и 
lin, ыр шге Шш; 
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Let us now divide both sides of (10) by 10% (e, x) and let x tend to infinity. Hence, 
by means of (7) and (13) we get 
1 
we? (67 , 
lim 2 | кыы и Es ола (17) 
€ 


aco у | ҰШЫР (ала) Е 


Making use of (16) and taking into account that the left-hand side of (17) does not 
vanish, it follows that 


с wt 2) (e, x) 502 (e, x) 
lim : 2 = li СЕБ И 
ае E 08) 


We note that relations (16) and (17), taken together, imply more than relations (16) 
and (18) taken together; from the former, we obtain some information about the 
type of singularity of the integrand in (17) (Dirac 6 type singularity). 

The relations (16) and (18) may be extended to W%” (e, x), namely 


IWEK (e, x) 
дє 


wi) (e, x) = 


whence 
1 
wD (д, з) 
THe" x) d 
DUAE WDE, a) 
2, € x 4” 
па ———————- — lim -цң-------- 19 
an WD (e, x) iss Y 19) 
we 1) (n, x) 
wt (е, x) 


where е! > =. From (7) and (16), it is seen that the numerator on the right-hand 
side of (19) vanishes and the denominator tends to infinity. Consequently 


WED (e', x) _ li SNS ШЕ 0 


Шал (а усы Due a) т (20) 
In similar manner we obtain 
(,2) (i, 2) 
in) ЕН EN аи 


gto PED (e, x) кю 90 (6, x) > 
From (6) and (21) follows immediately, that 


WEN (en) 2 SUD (e',x) 2 SHV Qe, x) 1 5%5(,4) 
EE de dum. eee n OD 0 (22 
ERU uo еп. 2) 
On account of the inequality (3), it is easy to obtain from (20), (21), and (22) 
the following formulae for P (e, x) 


P(e’, x) | 
EOS gu nba 23 
D PU on (23) 
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PC 2) (e, x) 
= 24 
a FEES T 
(4, 1) 
p DEEP I. 0 (25) 


lim соу ү = 
soo. P9 (e, x) 


ОЕ В. 
A discussion of results (20) and (21) is given in the last section of this paper 


3. Asymptotic formulae for the distribution function 


Besides P (e, x), defined in paper I and II, we introduce P(e, x) the probability 
of finding at least one particle of energy larger than ғ E, at depth x due to a primary 


particle of type i and energy Ев. 
We notice that we are able to define P^? as well as P in terms of the generating 
function С (e, U, us, x) (see Introduction of this paper), and so 


PED (e eee GU. ш a е 
=] (e, uy = б, иә = бі,, 5) (26) 
Р® (g, x) — 1 Фф) (e, О, О, x) 
In a manner similar to the case of РФР (в, x) (see paper I), we may show that 
PG (e, x) = eax VO (e, x) (27) 


where И (e, x) behaves іп a manner similar to ИФ 9 (e, x) (see relation (4)). Conse- 
` quently, using (2), we may at once write down the integral equations for PÓ (e, x): 


P (e, x) = SG (г, x) + se (ea) APOE = B, P® (e, x) (28) 


where 


А; Peau [ај ШЕ ЛЕС x — Ё) PO) (ео) 


(29) 


В; РФ (e, x) - fuf ш) (= ) 2 (7, x— £) PO? Е 3 s) Ра) Е 2 2 dé 
€ "~e 


Let us divide both sides of (28) Р® (e, x); we get 


m SEY (e, x) + SÈD (в, x) 


к — Ai—B (30) 


We shall now ow that 


lim A, = lim B; = 0 (31) 


aso å x-00 
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To prove (31) we write (29) in the form 
; wi 1) (n, x— 6) (2 | € | € (2) 
( m А ) ы (1) 15 ON 
24. cne f VG (e, x) 4 25 yV E € dé 
б 
„Гиз ,x—€) Е ІШЕ SPA 
PED E Io o y & s) yo کا‎ 3 е-е ЧЕ (39) 
0 
where we have used the notation (27) and (5a) and where 
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(7 Й Е 7 
E (33) 
D (s, е) ZE — (2) 
7 
0 0 
Because of the properties of w(^^ (e, x) we may write 
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soc UTE. M V vot y pi € dé 
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xoo0 $ 
lim | ИФ E г) yo - = s) e- € dE < со (35b) 
x = 0 8 xy 
0 


We shall show only that (35a) is true, since (35b) may be proved in quite a similar 
manner. From (35a) we have 
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in regard to properties of V“ (e, x). This proves (35a). Now we shall show that 


иф (ъа) 
36 
lin, 79 (s, x) = (36a) 


whe ® (qux). 
с ENAS 6b 
imie a T = (О) 

To do this it suffices to show that 
We 1) (%, x) E 
mE ОЕ 0 үе (37а) 


апа 


уэ) 
а 
since 10% 9 (e, x) is always non-negative. 
From (29) we see that for all 7, used there, we have 7 > 2e, thus, of course, 77 > €. 
It is easily to show that 


=0 ные (37Ь) 
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Thus, оп account of (38), we have | | 


S6D(nx) SHO (ny, x) 


< FI, PED,» | ue 


and ET 
| S62 (n. x) Se i E 
PO (e, x) р 
апа (6) we have — 
) г» = x ik 


КЕЗ 


Futher, on the dno of EE 


| 
t 


MM 


On the asymptotic Behaviour of Electron-Photon Cascade 


Let us write (41) in the form 


. ЭФ (а, 2) 9% 0 (e, x) 2 
PIE 4) | * $83 (e, x)f ^ 
= Then, because of (21), we have finally 


5% 2) (e, x) 
лз ACT A 
и P (e, x) 
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(42) 


Formulae (41) and (42) taken together enable us to obtain the үрне formulae 


for the distribution function Р® (п, mle, x) (see paper ID. 
We have 


_ PO (0,0, x) ~1 
Р® Olle, х) ~ 5%? (e, x) 
Р® (n, mle, x) = о (P (0,1, x) 


for pairs of numbers (n, m), different from (0,1), where n + m > 0. 


(43) 


In a very similar manner we obtain from (2), with help of (26), the integral equa- 


tions for P 9 (e, x) 


— 56 (e, x) 
= pàn (s, x) 
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The proof of (46) is very similar to that given above for (31). The only thing which 
needs some explanation is the proof of the relation 


(т, k) 
lim Б (Ёё 


x00 Vo 1) (е, x) _ 


for n > в and i, k, l, = 1,2. 
From (6), (20) and (22), using a method similar to that used in case of £ (36), we obtain, 


for Pe e 
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Thus, since (46) holds we DER finally from (44) 
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tends to infinity if the depth of the absorber increases to infinity independently of the 
type of primary particle initiating the cascade, provided the threshold energy of the 
apparatus registrating electrons and photons is the same. In other words, if the depth 
of the absorber is sufficiently large, the number of registrated photons is greater 
than the number of registrated electrons and this ratio will increase with increasing 
depth. Now relation (20) will „weaken“ the „strong“ relation (21). We interpret (21) 
as follows: | 

The ratio of the average number of photons of energy larger than &' E, to the 
average number of electrons of energy larger than е Ey, where е” > e, tends to zero if 
the depth of absorber tends to infinity independently of the type of primary particle. 
In other words: If the depth of the absorber is sufficiently large and the threshold 
energy for the photons is higher than for the electrons, then the number of registrated 
electrons is greater than the number of registrated photons and this ratio increases 
with increasing depth. 

As was said in the Introduction, our computations are done in the approximation 
A, i. e. for threshold energies chosen so that neither Greisen's photon effect (Greisen 
1949) nor the ionization has any significance. It is well known that these two factors 
act in a direction to stregthen the influence of photons and suppress the influence 
of the electrons (both, however, of low energy) at great depth, i. e., their action is in 


agreement with (21). But our result — as was already said — is independent of the 
two effects mentioned above. Therefore it seems to us that it would be interesting to 
verify our results experimentally for particles of higher energy, where Greisen's 
photon effect and the ionization do not obliterate our effect. We are, however, aware 
of some experimental difficulties. Although (21) holds, nevertheless for large but 
finite, depths, we do not expect the effect to be strong, because а) baryly differs 
from a(?, One may, however, compare the ratio of the numbers of photons and 
electrons for various depths. A further difficulty is that the number of photons as 
well as the number of electrons, falls off very rapidly with increasing energy for 
large depths. As already mentioned, this steepness of the average-number curve, 
considerably weakens the effect implied by (21). Consequently, one may expect 
that experiment would indicate that the effects of (20) and (21) blend into each other. 

Our results, as well as the results of paper II, are in agreement with the numerical 
data and curves for the ratio of first moments of the photon-electron cascade given by 
Jánossy and Messel (1951) and obtained by means of the saddle — point method. 
The derivatives of the functions represented by these curves are always positive. 
These function behave like eaVx where а > 0. 


The second main result is that 
P (0,0]e, x) ~ 1 
P (0,1|e, x) ~ 5% 2 (e, x) (43) 
PO (n, mle, x) = o (РФ (0,Це, а), 
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where п + т > 0 апа (п, т) + (0,1), or, written in terms of conditional probalities 
P® (0,1] e, x) 


Y: 3 PO (n, mle, x) 


п+т> 0 
РФ (n, m | e, х) 


23 > РФ (n, mle, x) E 


n+m>0 


Бо (01|, х) = 


--1 


(43 a) 
P (n, т | e, x) = 


We interpret (43a) as follows: 

if the depth of the absorber is sufficiently large and if we have registrated at 

least one particle, we may be sure that only one particle was registrated and 

that this particle was a photon. 
This statement means that the predominant part of registration events form the single 
particle registrations; other registrations may be neglected at large depth. 

In our considerations no attention was paid to the kind of the primary particle, 
initiating the cascade. In paper II, however, the role played by the primary particle 
was investigated exhaustively. Considering our résults in light of the results contained 
in paper II we are able to conclude as follows: 

if the depth of the absorber is sufficiently large and if we encounter at least 

one particle, we may be sure that there is only one particle. This particle is a 

photon belonging to a cascade due to a primary photon. 

Formulae (43) give the asymptotic solution of the distribution function for x tending 
to infinity. Therefore (43) is the first approximation for the distribution function for 
very large depth. In order to find the distribution function for a thinner absorber one is 
expected to find further approximations. These approximations are probably built 
up from higher factorial moments of the distribution. 

Eq. (48) may be interpreted in (41) and (42), which lead to (43). Thus: if at large 
depth, we have registrated at least one electron (photon) and an arbitrary number of 
photons (electrons), we may be sure that we have registrated only one electron (photon) 
and an arbitrary number. of photons (electrons). 

Finally we emphasize that in course of our work we did not make any assump- 
tions about the shape of the cross-section functions, except for the few general assum- 
tions noted in the Introduction. 

Our results are independent of the manner in which the particles in the cascade 
multiply. The only essential feature is that а) < a, 
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KPATKOE СОДЕРЖАНИЕ 


Конвент и Лопушанский, Несколько замечаний о асимптоматических свойствах 


‚ злектроно- фотоновой каскады для больших глубин абсорбера 


В работе показано, что: 
2 562 (е, а) Р 
п аали ДЛЯ в >E 
x00 SG, 1) (є, x) 
a также: 
50,2) (г, x) 
т — = 
SES SG, 1) (e, x) 


где: SG k) (e, x) обозначают первые статистические моменты для электроно-фото- 
нного каскада, притом, что функция распределения вероятности 


P (0,0 | e x) ^ 1 

PO (0,1 | e, x) ~ 8% 2 (e, x) 

PG (n, m | в, x) = o (PO (0,1 | e, x) 
(n, т) Æ (0,1); n + т > 0. 
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